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Abstract: The future of ground based gamma ray astronomy lies in large arrays of Imaging Atmospheric 
Cherenkov Telescopes (IACT) with better capabilities: lower energy threshold, higher sensitivity, better 
resolution and background rejection. Currently, designs for the next generation of IACT arrays are being 
explored by various groups. We have studied possible configurations with a large number of telescopes 
of various sizes. Here, we present the precision of source, shower core and energy reconstruction for 
gamma rays in the GeV-TeV range for different altitudes of observation. These results were obtained 
through tools that we have developed in order to simulate any type of IACT configuration and evaluate 
its performance. 



Introduction 

The future Imaging Atmospheric Cherenkov Tele- 
scopes (IACT) will be expected to discover new 
sources, enable more precise observations of 
known sources as well as contribute towards an- 
swering questions in adjacent fields like cosmol- 
ogy and particle physics. This will require large 
arrays of telescopes with lower energy thresh- 
old, higher sensitivity, better resolution and back- 
ground rejection. Here we present the performance 
of two possible array configurations for 7 photons 
in the GeV-TeV range at two different altitudes. 
This study was carried out through tools that we 
have developed and whose description is given in 
the next couple of sections. 

Simulation tools for the study of tele- 
scope designs 

IACT systems have a large number of parameters 
that can be optimised to improve detection capa- 
bilites, such as the number of telescopes, their po- 
sition, size and field of view and the altitude of ob- 
servation. In order to study IACT systems, we have 
developed a tool capable of simulating any type 
of telescope configuration and evaluating its per- 



formance. The atmospheric showers are simulated 
through the CORSIKA 1 package [1]. The reflec- 
tion of the Cherenkov photons from the shower by 
a parabolic mirror and their impact on the telescope 
camera are then carried out by our IACT simula- 
tion tool. The program allows complete freedom 
in the choice of telescope diameter, focal length, 
camera size, photomultiplier size, telescope posi- 
tion and orientation and altitude of observation. Up 
to 100 telescopes with variable individual charac- 
teristics can be simulated at the same time so as to 
enable the study of very large arrays. 

Shower parameter reconstruction 
methods 

The images obtained from the simulation tools 
can be used to reconstruct the source position, the 
shower's core position on the ground (herafter sim- 
ply called shower core or core) and its energy. 

Source reconstruction - In the camera frame of 
reference, each telescope gives a shower image 
whose axis contains the position of the source and 
if the images from several telescopes are super- 
posed (see fig. 1), then this source position corre- 
sponds to the intersection of the axes of all images. 

1. CORSIKA version 6.020 has been used here. 
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Energy reconstruction - The energy reconstruc- 
tion makes use of the linear relationship between 
the average number of photo-electrons obtained 
on a given telescope at a given distance from the 
shower core. Tables containing the average num- 
ber of photo-electrons for a wide range of distances 
and shower energies have been produced and are 
used to deduce the energy from the photo-electron 



content of telescope images. 



Figure 1: The superposed images of a 500 GeV 
shower obtained by four telescopes. The recon- 
structed axis of each image is shown in black. 



The source position can then be reconstructed by 
maximising the likelihood function: 
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Here is the content of the i th pixel from the 
j th telescope and dij is its distance from the im- 
age axis, (xij and y^) and (x c j,y c j) are the co- 
ordinates of the pixel and the centroid of the im- 
age, respectively. The following assumptions are 
made: (1) Each image axis is a straight line pass- 
ing through a point (x ,y ) common to all axes 
which gives the position of the source image in the 
camera frame of reference and whose coordinates 
are free parameters. (2) The distance of the pixels 
in an image from the corresponding axis (in other 
words the transvere image profile) is assumed to 
follow a Gaussian probability density function, at 
is the average standard deviation obtained from the 
Gaussian fit of transverse image profiles. (3) Each 
axis is made to pass through the centroid of the cor- 
responding image. The likelihood function is then 
maximised for x and y a yields the position of the 
source and reconstructed image axes as shown in 
fig. 1. 

Shower core reconstruction - In the ground frame 
of reference, the point of intersection of the axes 
from all shower images corresponds to the core 
position 2 . The core position can then be recon- 
structed through a likelihood minimisation similar 
to the one performed for source reconstruction. 



Array design 

Energy range and array parameters 

The physics goals and issues concerning 7-ray ob- 
servations depend on the energy domain being con- 
sidered 3 . The performance of current-day IACT 
indicates the following two energy domains 4 . 

High energy domain: 300 GeV - 10 TeV - 

Present-day telescopes have shown that this is the 
domain where IACT best operate and that good an- 
gular and energy resolution can be achieved using 
medium sized (10-15 m diameter) telescopes. The 
main objective in this domain is to further increase 
the sensitivity of the future telescopes . This can be 
achieved by spreading a large number of telescopes 
over a vast area. 

Low energy domain: < 30 GeV - As the shower 
size decreases with energy, the images have fewer 
photo-electrons and are more subject to fluctua- 
tions. This makes it harder to reconstruct shower 
parameters and separate 7-photons from hadrons. 
The main requirement in this domain is therefore 
the collection of a maximum amount of light from 
showers through the use of large sized telescopes. 
At the same time, since 7-ray fluxes tend to in- 
crease at lower energies, meaningful results can be 
obtained with only a few telescopes. 

We have chosen to work with a large number of 
12.5 m diameter telescopes aimed at the high en- 
ergy domain and 4 to 5 telescopes with large di- 



2. Strictly true only when the source is at the zenith. 

3. A discussion on the different energy domains and 
the scientific objectives can for instance be found in [2]. 

4. We have restricted this study to energies below a 
few tens of TeV. Beyond this limit, the flux from the 
sources become extremely low and very large surfaces 
of detection are required. 
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ameters of 30 m for the low energy domain. All 
telescopes have a moderate field of view i. e. 5.4°. 
The square shaped pixels on the camera have a side 
of 0.1°. The study is carried out at two altitudes 
above sea level : 1800 m and 3600 m. 

Optimum telescope separation 

Once choices concerning telescope size and num- 
ber have been made as a function of the energy do- 
main being targeted, the optimum inter-telescope 
distance can be determined. In order to do so, 
we have chosen to study the 7-ray response of a 
square unit of four telescopes by uniformly gener- 
ating vertical 7-rays over a large surface area. A 
simple trigger requiring that at least two telescopes 
have images with at least 50 photo-electrons is ap- 
plied and shower parameters are reconstructed for 
the passing events. 

This is carried out at 300 GeV with 12.5 m tele- 
scopes and 50 GeV with 30 m telescopes. These 
energies correspond to the lower ends of the high 
and low energy domain respectively 5 . 

The performance of the system for 7-rays, namely 
the efficiency of source, core and energy recon- 
struction is evaluated for different inter-telescope 
distances between 25 and 600 m. We specially fo- 
cus on the angular resolution shown here in fig. 2 
(top) as a function of the inter-telescope distance. 
The source position is best reconstructed when 
the inter-telescope distance lies within an opti- 
mum range of around 100-200 m. At higher alti- 
tude, this optimum range seems to be slightly nar- 
rower and the reconstruction of the source slightly 
poorer. The same optimum range can be found 
when looking at shower core reconstruction results 
(not shown here). The effective area (fig. 2, bot- 
tom) is essentially flat over a large range (100- 
300 m) of telescope separations with an optimum 
lying around 200 m. In both plots, one notes that 
the optimum range is independent of the energy 
range being considered and has only a slight de- 
pendence on the altitude of observation. 

Note that in order to focus on the intrinsic char- 
acteristics of the showers and avoid the introduc- 
tion of too many parameters we have neither added 
night sky background to the simulations nor per- 
formed image cleaning. The impact of the geo- 
magnetic field has been ignored as well. 
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Figure 2: The angular resolution of the four tele- 
scope unit (top) and its effective area (bottom) as a 
function of telescope seperation. 

Possible IACT array configuration 

The results shown in the previous section can then 
be used to obtain an array configuration. 

Low altitude configurations - The two configu- 
rations for low altitude observations are shown in 
fig. 3. In the first configuration, four telescopes of 
30 m diameter are placed at the corners of a 200 m 
sided square for low energy observations. This 
inter-telescope distance corresponds to the upper 
edge of the optimum telescope separation range. 
This choice allows us to optimise parameter recon- 
struction while keeping the largest possible effec- 
tive area, with only four telescopes. For the high 
energy domain a total of 33 medium sized tele- 
scopes are spread over a region with a radius of 
about 400 m. They are distributed in such a way 
that the inter-telescope distance in the resulting ar- 



5. One notes that while the low energy domain does 
not have a fixed lower limit, the combination of the small 
size of electromagnetic showers and the higher cosmic 
ray background levels for lower energies seem to in- 
dicate that gamma observation below 50 GeV will be 
highly problematic if not impossible. We have therefore 
chosen to carry out the inter-telescope distance optimi- 
sation at a relatively 'safe' energy, i. e. 50 GeV. 
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Figure 3: Array configuration 1 (left) and configu- 
ration 2 (right) at 1800 ma. s. 1.. 



ray is of 140 m. This lies in the middle of the 
optimum telescope separation range. In order to 
study the effect of a denser array, 16 more medium 
sized telescopes are added as shown in the right 
figure and the central telescope is replaced by a 
large sized telescope. The resulting configuration 
has an inter-telescope distance of 100 m which cor- 
responds to lower edge of the optimum telescope 
separation range. 

High altitude configurations - The two configu- 
rations are re-scaled 6 at higher altitudes so that the 
inter-telescope distance is 175 m for the large tele- 
scopes, 120 m for the medium sized telescopes and 
87 m for the denser configuration. 

Performance of possible arrays for 7- 
photons 

The performance of the two telescope configu- 
rations is studied by uniformly generating 7-ray 
showers over a surface of 2400m x 2400m at fixed 
energies and applying the simple trigger described 
earlier. Figure 3 shows the effective area and angu- 
lar and energy resolutions obtained for the two ar- 
rays. As can be expected, the effective area of the 
arrays at high altitude is smaller since the system 
has been re-scaled and the size of the Cherenkov 
ring is smaller. This effect is specially visible for 
higher energies where showers tend to get cut-off 
by the ground before having fully developed in the 
atmosphere. An angular resolution of around 0.07° 
is achieved at 1000 GeV. Note that while a four 
telescope system yields similar angular resolutions 
for shower falling within a radius of around 150 m, 
this angular resolution is calculated for all showers 
generated within a square region of 800m x 800m. 
Similar remarks can be made about the energy res- 
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Figure 4: The effective area and angular and en- 
ergy resolutions of the two IACT configurations at 
1800 m and 3600 m a. s. 1.. 



olution of around 7% achieved at 1 TeV. One also 
notes, that the use of a denser array (configuration 
2) does not seem to have any impact on the recon- 
struction of these parameters. Finally, the recon- 
struction capabilites of the arrays seem to improve 
slightly at lower altitude. Further details of this 
study can be found in [3]. 
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6. The system is scaled by a factor corresponding to 
the ratio of the Cherenkov ring size on the ground at 
both altitudes. The plots in fig. 2 suggest that there may 
be ways other than the simple re-scaling of the system to 
obtain a configuration at high altitude. However, the use 
of the same configuration allows us to have first compar- 
ison between systems at both altitudes. 



